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The electronic structure of the (2\/2 X \/2)R45 O/Cu(001) system has been cal- 
culated using locally self-consistent, real space multiple scattering technique based 
on first principles. Oxygen atoms are found to perturb differentially the long-range 
Madelung potentials, and hence the local electronic subbands at neighboring Cu 
sites. As a result the hybridization of the oxygen electronic states with those of its 
neighbors leads to bonding of varying ionic and covalent mix. Comparison of results 
with those for the c(2 x 2) overlayer shows that the perturbation is much stronger 
and the Coulomb lattice energy much higher for it than for the (2\/2 x y/2)RA5° 
phase. The main driving force for the 0.5ML oxygen surface structure formation on 
Cu(OOl) is thus the long-range Coulomb interaction which also controls the charge 
transfer and chemical binding in the system. 

PACS numbers: 68.43, 73.20.Hb, 71.15.Ap 

I. INTRODUCTION 

Experimental structural studies of Cu(OOl) with sub-monolayer coverage of oxygen shows 
clear signs of phase instability of the system (see for example Ref. p. The recent studies 
are a continuation of the debate that has circled around two main phases: a c(2 x 2) struc- 
ture on an unreconstructed Cu(001)i'§&! and a (2^/2 x \/2)-R45 structure which involves 
a reconstructed surfaced Some experiments have indicated the presence of both phasesSi. 
Still other have affirmed the presence of the c(2 x 2) overlayer at low coverages, but only in 
nanometer size domains and not in large well-ordered areasS Based on these experimental 
results, it is now widely accepted that for coverage lower than 0.34ML, oxygen forms c(2 x 2) 
islands. If the coverage exceeds this value, the (2a/2 x y/2)R45° structure is formedi. An 
array of experimental measurements based on surface X-ray diffraction!^, low energy elec- 
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tron diffraction0'0, surface extended X-ray adsorption!, and photoelectron diffractionlll have 
provided useful, and at times conflicting, information on (2v^2 x y / 2)i?45° superstructure of 
O on Cu(OOl). One of the common features of the superstructure, schematically illustrated 
in Fig. 1, from these studies is the missing of every fourth (010) row of atoms from the 
top Cu layer. Secondly, all studies point to the oxygen atoms occupying sites located on 
either side of the missing row between Cu atoms (denoted as Cu2 in Fig. 1). The structure 
undergoes further reconstruction (some displacement of Cu2 along the (100) direction) and 
relaxation of mostly the top Cu layer. There is, however, some discrepancy in the exact 
nature of the relaxation obtained from different experiments. For example, according to 
Ref. the heights zo, Zc u \ and zcui in Fig. 1, appear in the order zo > zq U 2 > Zcui, 
whereas Ref. [I] gives zq u \ > zo > zq U 2 (see Table | for details). 

The competition between the two observed phases of O on Cu(001) and the subsequent 
reconstruction of the surface with the (2-\/2 x y/2)R45 overlayer naturally beg the question 
about the origin of this complex behavior. The driving forces for such structural transfor- 
mation are expected to be related to the chemical binding and electronic structure of the 
system. Earlier theoretical studies, based on simple models, have correlated the structural 
properties of O/Cu(001) system to a strong pO - dCu hybridization^ and charge transfer^ 
in the vicinity of the surface. The rationale of these assumptions is supported by calcula- 
tions, which indicate the presence of the charge transfer from the metal atoms to oxygenl^E 
and a strong pO - dCu hybridizationEl for O/Cu(001) c(2 x 2) system. Our recent first 
principle calculations for the c(2 x 2) O overlayer on Cu(001)@ also confirm the importance 
of considerations of hybridization and charge transfer in the system and show that these fea- 
tures are eventually traceable to the long-range Coulomb interaction in the system. Clearly, 
the electronic structure and chemical binding in this system are complex and an adequate 
description of the relationship between the electronic and atomic structures of O/Cu(001) 
is yet to be established for the (2a/2 x y / 2)i?45° structure. Our goal here is to carry out 
such a study to understand the driving force for the observed surface reconstruction. For 
this purpose we take as benchmarks the two proposed surface configurations from Ref. [I] 



and [12] and compare their implications. As we shall see, with minor differences in the elec- 
tronic structure, both geometries provide the same rationale for the ensuing missing row 
reconstruction. Comparison of the results with those for the c(2 x 2) phase provides further 
evidence for the relative stability of the (2v^2 x y/2)R45° O overlayer for 0.5 ML coverage. 
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The calculations in the present paper have been carried out using the real space multiple 
scattering technique. We have chosen this particular method because apart from being 
based on first principles, it allows calculations of the charge transfer in the system with ease 
and provides a good estimate of the nature of the local bonding. This technique also lends 
itself nicely to the calculation of local densities of electronic states, their projections on the 
cubic harmonics, and the valence electronic charge densities. We also calculate the lattice 
Coulomb energies and the perturbation of the Madelung potential induced by the overlayer. 



II. COMPUTATIONAL DETAILS 



Our calculations are based on the density functional theory within the local density 
approximation (LDA)0 and multiple scattering theory. We use the local self-consistent 
multiple scattering (LSMS) scheme^, which is designed for extended heterogeneous solids. 
Within LSMS a compound is divided into overlapping clusters called local interaction zones 
(LIZs) centered around atoms of different local environments. For each LIZ, we solve a 
system of equations for the T— scattering matrix in the lattice site - angular momentum 
representation for the muffin-tin (MT) potentiaS 

E KKt - 4 n L \}T2i2 = 62%, (i) 

n\L\ 

where n denotes the lattice site number, L = l,m are angular momentums, t n i are the 
single-site scattering matrix elements, g^Li denote the free-electron Green's function matrix 
elements, and T^£,' are scattering path matrix elements. The solutions are used to determine 
the LIZ Green's functions, the local densities N^{E) of electronic states, their projections 
Nj^E) on the cubic harmonics, and local valence charge densities for the atoms located 
at the center of each LIZ. These charge densities are used further to solve the Poisson's 
equation for the entire extended system and hence to build a new potential for the next 
iteration of the self-consistent process. An application of this method to materials related 
to the work here has demonstrated its high efficiency and reliabilitylll'll'il. 

To take into account surface effects on the MT-potential, we follow the approach devel- 



oped in Ref. |24j. In this approach the space in the vicinity of the surface is divided into 
layers belonging to different atomic planes parallel to the surface and the interstitial charge 
density is supposed to be layer- dependent reflecting the asymmetry of the system. The 
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detailed description of the technique is given elsewhere^. Here we only remind our readers 
that in our calculations the Madelung potential includes monopole M^ J/3 , dipole and 
interstitial Vij[pj] termsll: 

y-l ad = E (ojpMZf, + d jP M™ jp + VM) , (2) 

where i,j denote the layer number, while a,/3 number the atoms in the 2D unit cell, djp is 
a dipole moment of the MT charge density, and the effective charge 

9jj8 = z op - qfp T + -^-pj R %- ( 3 ) 

The exchange and correlation parts of the potential are determined within LDA using 



the technique described in Ref. 25 



In this work the above formalism is applied to the case of (2^2 x V / 2) J R45° O overlayer 
on Cu(OOl). Fig. 1 is a schematic representation of the system. The parameters for the 
two proposed configurations from Ref. [I] and 12|, hereafter referred as "structure I" and 



"structure II", are summarized in Table |. The second and lower Cu layers are assumed to 
occupy bulk-like positions. The presence of oxygen atoms and the missing rows reduce the 
symmetry of the system and create nonequivalent environments for Cu atoms belonging to 
the same layer. We take this into account by allowing atoms of the first three layers to be 
nonequivalent as needed. They are marked by different numbers in Fig. 1. Thus we build 
LIZ's around oxygen and nine nonequivalent Cu atoms within the top five layers. We take 
the sites of missing Cu to be occupied by vacancies represented by MT-spheres with zero 
core electron density and requiring an additional LIZ to be built around the vacancy. These 
LIZ's contained 62 to 91 atoms depending on the local configuration, and the requirement 
that the calculated characteristics of the bulk system closely match those obtained from 
other reliable methods in the literature. 



III. RESULTS AND DISCUSSION 

Adsorption of gases on surfaces causes both short- and long-range effects on the electronic 
structure of the system. The short-range effect originates from changes in the local envi- 
ronment of some surface atoms which lead to modifications of their chemical bonds. The 
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long-range effect induced by the Coulomb interaction between the chemisorbed overlayer 
an the substrate, on the other hand, may cause differing shifts in the local potentials at 
atomic sites with differing environments. As a result, charge transfer may be induced and 
the electronic state hybridizations at these sites may be modified. 

Here we examine in details both the short- and long-range effects of the chemisorption of 
an oxygen overlayer on Cu(OOl) with the view of understanding the driving forces responsible 
for the observed missing row reconstructed (2\/2 x \/2)i?45 superstructure. We present first 
a comparative study of the short-range effect in structures I and II through examination of 
the local densities of electronic states at the oxygen and top Cu sites and their hybridizations. 
To gain further insight into the nature of the O-Cu bonding and subtle differences emerging 
from structures I and II, we resort to an analysis of the projections of these densities of 
states on cubic harmonics. These short-range effects, however, are related to the long-range 
ones responsible for shifts in the local potentials at each site. We do this next through a 
systematic evaluation of the effects of the Coulomb interaction on local potentials. This is 
followed by an analysis of the charge transfer in the system. Finally, we comment on the 
relative stability of the c(2 x 2) phase, in view of the findings here for the (2\/2 x y / 2)i?45° 
phase. 



A. Local densities of the pO- and dCu-electronic states 

The calculated densities of the p-states of the oxygen and d-states of its nearest neighbors 
Cul, Cu2 and Cu3 atoms for structures I and II are plotted in Fig. 2. The densities of the 
dCuA through <iCu9 states (N^ U (E)) are shown in Fig. 3. Since the difference between 
Nj u (E) of Cu6 and Cu7 is found to be negligible, only one is displayed. Figs. 2 and 3 
indicate a dramatic difference in N™(E) for atomic sites with different environments. The 
Nj u (E) of the Cu6 and Cu8 atoms belonging to the third and fourth layers, respectively, 
differ slightly from bulk values, while that of Cu9 coincides with that of the bulk. The <iCu4- 
and <iCu5-states, on the other hand, consist of a pronounced narrow peak together with some 
low energy structure. Finally, the densities of the pO-states and d-states of oxygen's nearest 
neighbors are found to be significantly split and structured. This splitting results in two 
substructures ("a" and "b" in Fig. 2), which have the same energetic position and similar 
shape for all spectra. Such a behavior clearly indicates a strong hybridization and covalent 
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binding between the ^O-states and the d-states of the neighboring Cu atoms. A relatively 
high intensity of the low energy peaks "a" in N^ u3 (E) suggests that the c/Cu3-states are 
more involved in the p — d hybridization than the dCul- and dCu2-states. To understand the 
origin of this difference we should remember that the degree of hybridization of electronic 
states depends on their spatial and energetic overlap. Since, in the present structures, the 
0-Cu3 interatomic distance is longer than O-Cul and 0-Cu2 one the strong pO-dCu.3 
hybridization does not result from the larger spatial overlap. The energetic overlap, on the 
other hand, can be quantitatively characterized within the multiple scattering theory by 
comparison of the single-site resonance energies (E™ es ) for neighboring sites, which reflect 
energetics of the local potentials and primarily govern the energy of local subbands@. The 
calculated E™ es for structures I and II, in Table |J, show substantial variation of E™ es from 
site to site. As we shall see in section C, these variations are a result of the long-range 
interaction and charge transfer in the system. Since from Table ||, the energetic separation 
between the pO and g?Cu3 resonances is smallest, we can conclude that the strong covalent 
pO - dCu3 binding is caused by a large energetic overlap of these states. Also, as -E^" 4 
and Egf are much higher than the others in Table 0, there is only a small coupling of the 
dCuA- and c?Cu5-states with those of their neighbors as indicated by the low intensity of the 
structures seen in the low energy region in Fig. 3. 

Having established the strong coupling between the pO and <iCu3 electronic states, we 
proceed next with an analysis of the projections of the local densities of states on cubic 
harmonics (N^^E)) which provide more detailed description of covalent binding in the 
system and help isolate the reasons for the a noticeable difference between the shape and 
splitting of for structures I and II. 

B. Projections of the densities of states on the cubic harmonics 

Since in the geometric structures under consideration, the local environment of the oxygen 
atoms have a reduced symmetry (see Fig. 4) because of the missing rows, the local electronic 
states are expected to be anisotropic. Therefore, to build a consistent picture of chemical 
binding in the system, it is important to analyze the projections of electronic state Nl^(E) 
on the cubic harmonics. Plots in Fig. 5 of such projections of pO- and the five <iCu3-states, 
calculated for structure I, indicate that they differ dramatically from each other. The O- 
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p z and Cu3-rf 3z 2_ 1 states are found to be strongly split. They form a common subband 
represented by peaks "a" and "b" , which are separated by about 0.5Ry. This strong splitting 
together with the involvement of major part of these states in hybridization (redistributed 
into " a" and " b" regions) are a clear signature of strong covalent binding. The 0-p x and Cu3- 
d xz as well as 0-p y and Cu3-d yz states also form common split subbands, but not as strongly 
as the 0-p z - Cu3-ds z 2_ 1 states. The cause of this difference in the hybridization becomes 
clear through a comparison of the angular distribution of corresponding cubic harmonics, 
which shows that the 0-p x - Cu3-d zx and 0-p y - C\x3-d zy spatial overlap is lower than that 
of 0-p z - Cu3-c?3 2 2_ 1 (Fig. 5). The C\x3-d xy and C\x3-d x 2_ y 2 states cannot be coupled with 
the oxygen electronic states by symmetry. Instead they hybridize weakly with the electronic 
states of Cu4 and Cu5, as reflected by low intensity structure at high energies in Fig. 5. 
The weakness of this particular hybridization can be traced to the high energetic separation 
of the single-site resonances (see Table H). 

In Fig. 6, the corresponding splitting of the 0-p z - Cu3-d 3z 2_ 1 subband for structure 
II is found to be larger (5.7 Ry versus 5.0 Ry) than that for structure I, implying stronger 
covalent binding. Although for the 0-p x - C\x3-d zx and 0-p y - Cu3-d zy subbands it is harder 
to quantify the splitting, Figs. 5 and 6 indicate that these subbands are also split stronger 
for structure II than for structure I. Again, these differences in splitting can be traced to 
the energetic separation between the pO and <iCu3 single-site resonances found for the two 
structures (see Table |I|). 

The projections N^{E) calculated for the other nearest neighbors of oxygen, namely, 
Cul and Cu2 are shown in Fig. 7 and 8, respectively. Only a minor part of the <iCul- 
and <iCu2-electronic states are found to contribute to the peaks responsible for the pO - 
dCu hybridization. As mentioned above, the single-site <iCul- and <iCu2-resonances are 
energetically separated from the pO-resonance that reduces the hybridization. The low 
energetic structures of N^(E) and N^ 2 {E) align with "a" peaks of Np m (E), which are 
formed with the pO-dCu3 hybridization. This suggests that the c/Cul- and dCu2-states just 
admix with the pO-<iCu3-subband and the O-Cul and 0-Cu2 covalent bonds are weaker 
than the 0-Cu3 bond. 
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C. Long-range Coulomb interaction and charge transfer 

For a deeper understanding of the factors controlling the electronic and atomic structures 
of the O/Cu(001) system, we turn now to considerations of long-range effect of the Coulomb 
interaction as expressed by the Madelung potential and the subsequent charge transfer in 
the surface region. To get a simple and intuitive picture we first calculate the Madelung 
potential for several surface structures of Cu(OOl) with 0.5 ML of O, within a simple "frozen" 
charge density approximation using Eq. 2. Namely, we calculate the effective charges q ia 
introduced in Eq. 3 using bulk-like electronic density for Cu and atom-like density for O. 
Such an approach gives us a pure geometric effect of the long-range interaction. As expected, 
the presence of the surface itself causes only small perturbation in the Madelung potential. 
This is represented by the low values in the column under "clean" (Cu(001) in Table |TJ. 
The admission of 0.5 ML oxygen on bulk terminated Cu(001) causes large perturbation in 
the Madelung potential at atomic sites in the top few layers, as seen in the column denoted 
by BT in Table |T|. By removing every fourth row of Cu atom from this bulk terminated 
O/Cu(001) system, we find the perturbation of the Madelung potential to be significantly 
reduced in the top Cu layer and slightly increased in some of the Cu sites below (see column 
labeled MR in Table [TTTD . The structures BT and MR are naturally artificial not only in 
terms of the positions of the top layer Cu atoms, but also that of the O which is taken to 
occupy the hollow sites in the surface plane. They have been chosen to provide benchmarks 
for the relative stability of the two observed O superstructures. The calculated perturbation 
of the Madelung potentials for structure II, i.e. (2^/2 x y / 2)i?45° O phase with Cu atoms in 
relaxed and reconstructed positions, in Table |T|, indeed display significant reductions from 



the values of BT and MR configurations, for all atomic sites. Finally, in Table ITJ we have 
also included the calculated quantities for the c(2 x 2) structure (parameters taken from 
Ref. |I] for a coverage less than 0.34 ML) for comparison. Clearly, the c(2 x 2) phase, albeit 
for 0.5 ML coverage, invokes a much stronger perturbation of the Madelung potential than 
the one in structure II. 

An electronic response to the strong perturbation of the Madelung potential is expected 
to be charge transfer in the vicinity of the surface. The resulting charge distribution would 
then determine the degree of ionicity of the system, and would influence single-site resonance 
energies, and consequently, the formation of covalent bonds. Charge transfer is thus a very 
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important characteristic of systems. To characterize it quantitatively, the system has to be 
divided into spaces belonging to each atom and the charge within this space calculated. In 
the framework of the MT- approximation, non-overlapping spheres centered on the atomic 
sites can naturally represent these spaces. A complexity of the present system is that 
interatomic bond lengths for a number of atoms are different because of the differences in 
their local environment. Thus different volumes belong to different atomic sites. To make 
the site charges comparable, we have built the spheres with the fixed radius equal to one-half 
of the shortest interatomic distance (1.807 A) and calculated the local valence charges by 
integrating the self-consistent valence electronic density over these spheres. The results in 



Table [TV] show a significant electronic charge transfer from Cu to O in both structures I 
and II. The charge transfer pattern is, however, interesting as Cu2, Cu4 and Cu5 donate 
part of their charge but Cul and Cu3 do not. These latter two nearest neighbors of O, in 
fact, gain a small amount. Such a complex character of the charge transfer reflects the local 
potential perturbations caused by the missing Cu rows induced by the oxygen overlayer. 
Indeed, from Table |TTTj and Table |Iy| it follows that sites loose electronic charge if the initial 
Madelung potential is higher and gain it if the potential is lower. Table [TV] also reflects 
subtle differences in the charge transfer to the different sites in the two structures (I and II) 
considered here. Recall that covalent bonds have also been found stronger in structure II as 
compared to structure I. 

Naturally the Coulomb interaction of the different sublattices plays a very important role 
in the above considerations. To get further insights into the difference arising from the two 
structural models chosen in this paper, we compare the Coulomb lattice energies E\ at for 
these O/Cu(001) phases. Within the MT approximation!! E[ at can be expressed as follows: 

Elat = ^ 53 M ia,jpqia<ljf3 + -^rO- 2 { U iaPi) 2 + ^QiaPi^ia) , (4) 
i,a i,a MT 

where r^ T and Ui a are the MT-radii and MT-sphere volumes, respectively. The first sum 
of Eq. 4 represents the Madelung energy and the second comes from the charge neutrality 
condition. The E[ at thus includes the Coulomb interaction between the MT charges and that 
between them and the layer-dependent uniform interstitial charge density in the system. For 
comparative purposes it is sufficient to include the contribution of the oxygen overlayer and 
only the top five Cu layers to Ei at . Eq. 4 then yields values of -13.864Ry and -15.420Ry for 
Ei at , per 2D unit cell, for structure I and structure II, respectively. The lower value of Ei a t 
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along with the finding of stronger O-Cu covalent binding in structure II than in structure I, 
indicate preference for the former structure. 

We have also calculated the Coulomb lattice energy for the c(2 x 2) O overlayer on Cu(001) 
using our self-consistent calculation results from Ref. and obtained the value of E\ at =- 
8.582Ry. This is much higher than that in the case of both structures I and II. This result 
is not surprising as we have already shown in Table |TTT| that the perturbation to the local 
Madelung potential at the atomic sites is also larger for the c(2 x 2) superstructure. Based 
on these findings we conclude that the long-range Coulomb interaction is the main factor 
which controls the phase formation in the O/Cu(001) system and makes the (2\/2x v^2)-R45° 
phase preferable at oxygen coverage of 0.5ML in agreement with experimental findings. 



IV. CONCLUSIONS 



The results presented above, provide the following microscopic picture for the formation 
of the (2y/2 x y / 2)i?45° O overlayer on Cu(001). Oxygen adsorption in the four-fold hollow 
sites of Cu(001) induces a strong Coulomb perturbation on the top layer Cu sites. Missing 
row reconstruction of the Cu surface, lattice relaxation, and charge transfer help signifi- 
cantly reduce this perturbation. In response to the Coulomb perturbation, the electronic 
charge transfer from the host surface to O is found to have a complex character, involving 
two of the four nearest neighbors and two next nearest neighbors of O. The perturbation is 
screened partially and the self-consistent local potentials are such that the energetic separa- 
tion between the single site pO- and c/Cu-resonances is different for different Cu sites. This 
energetic separation essentially determines the degree of the pO-dCu hybridization and the 
strengths of the covalent bonds formed between oxygen and its nearest neighbors. Thus, 
the Cu3 atom located right under O is involved in strong covalent binding with O, whereas 
weak covalent bonds are formed between the O and its other nearest neighbors. The O- 
Cu3 binding is found to be stronger in structure II than in structure I. Furthermore, as 
the Coulomb lattice energy is lower for the structure II than for structure I. In extending 
these consideration to the observed c(2 x 2) O overlayer on Cu(001) (for coverage less than 
0.34ML), we conclude that the c(2 x 2) structure is unstable, for half monolayer coverage, 
because it induces a strong perturbation at the surface that increases the Coulomb lattice 
energy of the system. Thus, according to our results, the long-range Coulomb interaction is 
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the main driving force for the O/Cu(001) surface atomic structure formation. 

Finally, we note that adsorption processes on metal surfaces are extensively studied by 
means of such methods as Monte Carlo or molecular dynamics simulations. Many of them are 
based on the model potentials, which depend only on the local atomic environment, namely 
on the interatomic distance and the number of the nearest neighbors. Our findings indicate 
that an application of such approaches to O/Cu(001) can hardly be successful. Indeed, 
according to our results, the long-range interaction in the system governs not only charge 
transfer (as usual), but also such "short-range" effects as the electronic state hybridization 
and covalent bond formation. Such an implicit impact cannot be taken into account even 
if a model potential somehow includes the long-range interaction along with a traditional 
short-range part. This suggests that a proper model potential for O/Cu(001) can be built 
either by parameterization of the energetic surface obtained from first principle calculations 
or by development a sophisticated form, which takes into account the interplay between the 
long- and short-range parts. 
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Figure Captions 

Fig. 1. Schematic illustration of the O/Cu(001) (2y/2 x y/^)RA5° structure. 

Fig. 2. Densities of the p-electronic states of oxygen and d- electronic states of its nearest 

neighbors calculated for structure I (solid line) and structure II (dashed line). 

Fig. 3. Densities of the ^-electronic states of the Cu4 through Cu9 atoms calculated for 

structure I (solid line) and structure II (dashed line). 

Fig. 4. Local environment of oxygen atom in the (2y/2 x y / 2)i?45° O/Cu(001) structure. 
Fig. 5. Projections of the pO- and a?Cu3-electronic states calculated for structure I. 
Fig. 6. Projections of the pO- and a?Cu3-electronic states calculated for structure II. 
Fig. 7. Projections of the rfCul-electronic states calculated for structure I (solid line) and 
the structure II (dashed line). 

Fig. 8. Projections of the <iCu2-electronic states calculated for structure I (solid line) and 
structure II (dashed line). 
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1.94 
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°the Cu2 atom displacement along 
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TABLE II: Energetic positions (Ry) of the single site pO- and <iCu-resonances E™ es counted from 


the MT-zero level 






Site 


Structure I 


Structure II 





0.232 


0.222 
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0.354 


0.377 
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0.416 


0.447 
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0.316 


0.283 
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TABLE III: The variation of the Madelung potential (in Ry) at different atomic sites with respect 



to the bulk value for 


the following configurations: clean Cu(OOl) surface (Clean) 


i; O on bulk 


terminated Cu(OOl) (BT)); (2y/2 


x y/2) R45° O 


on bulk terminated O/Cu(001) with 


missing rows 


(MR); structure II; c(2 


x 2) on 


Cu(001) (c(2 


x 2)) 
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-0.207 
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TABLE IV: 


Variations of the local valence electronic charge from the bulk-like charg 


e for Cu sites 


and from the atomic-like charge for the site. 
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